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ABSTRACT 

We consider the formation of double white dwarfs (DWDs) through dynamical interactions in globular clus- 
ters. Such interactions can give rise to eccentric DWDs, in contrast to the exclusively circular population 
expected to form in the Galactic disk. We show that for a 5-year Laser Interferometer Space Antenna (LISA) 
mission and distances as far as the Large Magellanic Cloud, multiple harmonics from eccentric DWDs can be 
detected at a signal-to-noise ratio higher than 8 for at least a handful of eccentric DWDs, given their formation 
rate and typical lifetimes estimated from current cluster simulations. Consequently the association of eccentric- 
ity with stellar-mass LISA sources does not uniquely involve neutron stars, as is usually assumed. Due to the 
difficulty of detecting (eccentric) DWDs with present and planned electromagnetic observatories, LISA could 
provide unique dynamical identifications of these systems in globular clusters. 

Subject headings: Stars: Binaries: Close, Stains: White Dwarfs, Gravitational Waves 



1. INTRODUCTION 

Double white dwarf binaries (DWDs) are the single most 
abundant and a guaranteed class of gravitational wave (GW) 
sources for the Laser Interferometer Space Antenna (LISA; 
Bender et al. 1998). Their formation encompasses two mass- 
transfer phases, with at least one leading to a common- 
envelope phase shrinking the orbit to periods of 0. 1-1000 hr 
(Han 1998; Nelemans et al. 2001). Once formed, orbital an- 
gular momentum losses due to GW emission and tidal inter- 
actions continue to shrink the orbit until the least massive WD 
fills its Roche lobe and transfers mass to its companion. De- 
pending on the binary mass ratio, the mass transfer leads to a 
fast merger or a long-lived phase of orbital expansion (Marsh 
et al. 2004; Gokhale et al. 2007, and references therein). 

As the progenitors of DWDs are tight enough to undergo 
multiple mass-transfer phases, tidal forces are usually thought 
to circularize their orbits long before they form a DWD. 
DWDs are therefore expected to be born with circular or- 
bits. However, dynamical interactions in globular clusters 
(GCs) form DWDs at enhanced rates compared to the field 
(Shara & Hurley 2002). Through these interactions, DWDs 
can form with non-zero eccentricities; these can subsequently 
be damped by GW emission, but can also be excited by fly-by 
encounters with single and binary stars. Contrary to common 
assumptions, the identification of eccentric stellar-mass bina- 
ries in the LISA GW data stream therefore does not uniquely 
identify them as neutron star (NS) binaries. 

Due to their inherent faintness, DWDs are difficult to de- 
tect at GC distances or beyond, even with the SDSS, HST, 
or JWST in the future. LISA could therefore provide unique 
dynamical identifications of eccentric DWDs in GCs. As the 
potential of studying WD interiors through tidal effects in- 
creases dramatically with increasing orbital eccentricity, ec- 
centric DWDs offer a unique opportunity to study stellar as- 
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trophysics with LISA. 

In this Letter, we explore the formation of DWDs with non- 
zero eccentricities (e > 0.01) and orbital periods P < 5000 s 
through dynamical interactions in GCs, and assess their de- 
tectability by LISA. Since mass-transferring binaries are not 
expected to retain appreciable eccentricities for long periods 
of time, we focus our attention on detached DWDs. While 
the importance of stellar-mass binaries in GCs for LISA has 
been discussed previously (Benacquista 1999, Benacquista et 
al. 2001, Benacquista 2001), the possible presence of eccen- 
tric DWDs in the LISA band has not been considered before. 

2. DOUBLE WHITE DWARFS IN GLOBULAR CLUSTERS 

Ivanova et al. (2006) recently explored the formation and 
evolution of WD binaries in GCs using an updated version 
of the Monte Carlo code described in Ivanova et al. (2005). 
In this study, we adopt these sample simulations to address 
the relevance of eccentric DWDs as LISA sources. The code 
combines a simple two-zone core-halo GC model with the 
STARTRACK binary population synthesis code (Belczynski et 
al. 2007) and the FEWBODY small N-body integrator (Fregeau 
at al. 2005) to treat dynamical interactions between single 
stars and binaries. 

In the simulations, the two most important types of stellar 
interactions creating eccentric DWDs in the LISA band are 
exchange interactions and physical collisions. Exchange in- 
teractions are a subset of the broader class of binary-single 
and binary-binary interactions. Eccentric LISA DWDs form 
through this channel by exchanges between single WDs and 
binaries with a WD component. During the exchange, the 
two WDs become bound, while the companion to the initial 
binary WD is ejected from the system. The resulting binaries 
have average post-exchange eccentricities e ~ 0.8. Alterna- 
tively, physical collisions between WDs and red giants form 
eccentric DWDs when the giant's envelope gets disrupted and 
its degenerate core becomes bound to the initially single WD 
(Lombardi et al. 2006). LISA DWDs formed this way have 
average post-coUosion eccentricities e ~ 0.5. We stress that 
the quoted average eccentricities are for LISA DWDs only 
(average post-exchange and post-collision eccentricities for 
the entire DWD population are lower). 

In addition to direct formation, we find that non-zero ec- 
centricities of DWDs can be significantly enhanced through 
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TABLE 1 

Globular cluster model parameters. 



TABLE 2 

Statistical PROPERTIES of eccentric LISA DWDs. 



GC model 




frh 


Z 






GC model 


{N) 


Pi 


Pi 


^3 


Typical 


105 


1.0 


0.005 


10 


40 


Typical 


0.39 (0.01) 


0.32 (0.01) 


0.06 


0.01 


47 Tuc-like 


2.5 X 10^ 


3.0 


0.0035 


11.5 


57 


47 Tuc-like 


1.9(0.03) 


0.84 (0.03) 


0.55 


0.29 


Terzan 5-like 


8 X 10 '' 


0.93 


0.013 


11.6 


49 


Terzan 5-like 


1.6(0.03) 


0.80 (0.03) 


0.48 


0.22 



Note. — GC core number density (pc ), half-mass relaxation time 
frh (Gyv), metallicity Z, one-dimensional velocity dispersion (T| (kms"'), and 
escape velocity Ve^c (kms-'). 

one or multiple binary-single star scattering events. These in- 
teractions are either distant fly-by encounters or close, strong, 
and possibly resonant interactions. Even for distant encoun- 
ters (at more than 5-10 times the binary separation) and ini- 
tial eccentricities e < 0.1, Heggie & Rasio (1996) find that the 
relative change in eccentricity can be as high as 30%. Con- 
sequently, the absolute eccentricity change in a given inter- 
action is insignificant for circular binaries, but increases with 
the binary's eccentricity and the cumulative effect of multiple 
interactions. On the other hand, close and possibly resonant 
interactions are rare, but can change the binary eccentricity 
drastically in just a single encounter 

The possibility of forming eccentric DWDs through the 
above channels is in stark contrast to DWDs in the Galactic 
disk which are expected to be invariably circular Despite the 
many dynamical interactions, circular binaries still dominate 
DWD populations in GCs, although we will show in the next 
section that the number of eccentric DWDs is non-ngligible. 
Circular DWD binaries could in principle also become eccen- 
tric via binary-single star encounters, although the efficiency 
of the process is small (Heggie & Rasio 1996). Nevertheless, 
all possible eccentricity inducing effects are included in the 
simulations described below, regardless of their relative effi- 
ciency. Since the detectability and properties of GW signals 
from eccentric DWDs are independent of the details of their 
formation, we do not distinguish between the different forma- 
tion channels in the remainder of this paper 

3. MONTE CARLO SIMULATIONS 

We use simulations by Ivanova et al. (2006) to estimate 
the expected number of eccentric DWDs in the LISA band 
(10~'*-10~' Hz). The outcome of the simulations depends on 
the adopted GC model parameters (e.g., cluster core density, 
metallicity) and stellar and binary evolution parameterizations 
(e.g., initial masses and orbital elements, efficiency of orbital 
angular momentum loss mechanisms). We consider three GC 
models representative of a typical Galactic, 47 Tuc-like, and 
Terzan 5-like GC (see Table [T]i, and GC ages in the range 
of 10-13Gyr. All stellar and binary evolutionary phases are 
treated as in the standard model of Ivanova et al. (2006). 

The maximum number of eccentric DWDs present in the 
LISA band at any given time in the 10-13 Gyr time inter- 
val is 2, 8, and 7, for the typical, 47 Tuc-like, and Terzan 
5-like GC models, respectively. Using the total number of 
systems formed in the 10-13 Gyr time window, we calculate 
the formation rates of eccentric LISA DWDs in the typical, 
47 Tuc-like, and Terzan 5-like GC models to be ~ lO'^^yr"', 
~ 3 X 10"** yr"', and ~ 5 x 10"^ yr"'. For each GC model, the 
probability that eccentric LISA DWDs are present at any 
given time in the 10-13 Gyr age range is furthermore well fit- 
ted by a Poisson distribution with mean (A^) listed in Table |2l 
It follows that the probability to find at least 1 eccentric DWD 
in the LISA band is 32% for the typical Galactic, 84% for the 



Note. — Mean number (N) of eccentric LISA DWDs and Poisson proba- 
bilities Pn that at least N such DWDs ai'e present at any given time. Numbers 
in parentheses correspond to DWDs with at least two harmonics detectable 
at a distance of 10 kpc with SNR > 8, for a 5 year LISA mission and a single 
time-delay interferometry observable. 

47 Tuc-like, and 80% for the Terzan 5-like GC model. 

The hatched probability distribution functions (PDFs) in 
Fig. [U show the statistical properties of the eccentric LISA 
DWDs in the three GC simulations. For all three models, the 
number of systems decreases monotonically with increasing 
orbital frequency v = 1/P. The abrupt drop in the number 
of systems at logi^ < -3.7 corresponds to the adopted cut-off 
orbital period of 5000 s for sources in the LISA band. The 
orbital eccentricity distributions (affected by post-formation 
GW emission) typically peak at e ~ 0.1 and decrease slowly 
toward higher values. In all three cases, there is a non- 
negligible fraction of systems with e up to ~ 0.7. The ma- 
jority of systems consists of at least one WD more mas- 
sive than ^ Q.6Mq. The total mass distributions therefore 
peak at relatively high masses Mi +M2 — 1.O-1.2M0. The 
chirp mass distributions peak at ~ O.3-O.4M0. The typ- 
ical merger lifetime Teyo of an eccentric LISA DWD is 10- 
100 Myr. However, a substantial number of systems also have 
considerably shorter lifetimes down to ~ 1 Myr. 

4. IMPLICATIONS FOR LISA 

Gravitational waves from binaries with non-zero eccentric- 
ity produce radiation at frequencies that are integer multiples 
of the orbital frequency. To compute the LISA sensitivity to 
eccentric DWDs, we model GWs at the leading quadrupole 
order The (angle-averaged) SNR at which a given binary can 
be detected over an observation time Tabs is then given by 



(s/n)>e((s/n);)=e/ 



where n labels the harmonics at frequency /„ ~ n z/. 



dlnf„, 



/ic,„ = (nd)-' {2E„/f„) 



1/2 



(1) 



(2) 



is the characteristic amplitude of the «-th harmonic, and 
'^?ms(/n) = i'^^/^)Sh{fn)fn is the root-mcan square value of the 
noise averaged over source position and orientation; Sh(f„) is 
the one-sided noise spectral density, which we computed us- 
ing Larson's online sensitivity curve generator''. In Eq. (|2]), 
En is the time derivative of the energy radiated in GWs at fre- 
quency /„ ~ nv, which to lowest order is 



£„ = (32/5)(27ri/M)"'/^g(«,e) 



(3) 



Here, g(n,e) is a function of the orbital eccentricity given by 
Eq. (20) in Peters & Mathews (1963). The term /„ ~ ndv/dt, 
to the leading quadrupole order, is 

f„ = 96{\QiT)-^ M^l^il-Kuy^l^ F{e) , (4) 

* http://www.srl.caltech.edu/~shane/sensitivity/ 
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— Statistical properties of eccentric DWDs (e > 0.01) in the LISA band (P < 5000 s). Hatched histograms represent the entire eccentric DWD 
grey shaded histograms represent the subpopulation of systems for which at least 2 harmonics are individually observable by LISA at SNR > 8, 
5 year LISA mission, a distance of lOkpc, and a single time-delay interferometry observable. As the PDFs are normalized to unity, the vertical scale 
panel to panel. 



where f (e) = [ 1 + (73 /24) + (37 /96) e^]/( l- e^^^ . The lat- 
ter function tends to unity for e ^ 0, but can become of the 
order of a few for e-values typical of the systems produced in 
our simulations, e.g. F(0.2) = 1 .29 and F(0.5) = 4.88. 

Using Eqs. ([TJ-lllll, we compute the number of detectable 
eccentric LISA sources, requiring that at least two harmonics 
be individually observable at SNR of 8 or higher. This thresh- 
old is conservative, since we still do not exactly know how the 
data analysis will be carried out (note that in the LISA litera- 
ture a SNR = 5 threshold is usually considered). The choice 
of assuming that at least two harmonics are observable sepa- 
rately, is dictated by reasons of robustness: from an analysis 
point of view, an eccentric DWD looks like a collection of n 
"circular" (one single harmonic) binaries, with suitable am- 
plitudes and frequencies. There is growing evidence that data 
analysis schemes that are already in hand for circular binaries 
(cf. e.g. Arnaud et al. 2007, and references therein) will per- 
form well at such moderate-to-high SNR. We can therefore be 
confident that even an analysis, comprising of a first pass of 
the data to search for "circular" binaries, followed by a second 
stage in which one looks for a more consistent fit to the data 
by allowing that some signals are actually harmonics of the 
same system is available. In practice, it is quite conceivable 
that a full-blown analysis allowing for non-negligible eccen- 
tricity can be carried out from the beginning. In this case the 
number of detectable systems per GC presented here is con- 
servative and should be regarded as a lower limit. 

In Fig. m we show the detectability results for the typical 
Galactic GC model, adopting a distance of lOkpc. The am- 
plitude on the vertical axis is the characteristic amplitude hc.„ 
multiplied by the square root of nm\[\,df„/dt x (Jobs/ f„y\ to 
consistently account for the frequency band swept by each 
source (all the systems in our simulations satisfy the condition 
dfn/dt X (Tobs/Zn) < 1 for Tobs = 5 yr). The height of each dot 
above the noise level then represents the relevant SNR value 
[cf. Eq. ([T]i]- We consider Tobs = 5 yrs, the current require- 
ment for the minimum mission lifetime. We also use only one 
Michelson observable in the analysis; if all six Doppler links 
are available over Tobs, the SNR would increase by a factor 
~ ^/l. Note that only systems for which at least two harmon- 
ics have ((S/N)^)'/^ > 8 are plotted in the figure. 

In Table [3] we list the number of detectable eccentric 
DWDs formed in the 10-13 Gyr age range, for distances 
d = 3, 10,50kpc. Depending on the GC model, the numbers 




frequency / Hz 



Fig. 2. — LISA sensitivity to eccentric DWDs in the typical Galactic GC 
model. The solid lines show the angle-averaged level of the rms noise (thin: 
instrumental noise; thick: total instrumental and galactic foreground noise). 
Squai'es indicate the amplitudes of DWDs for which at least two harmonics 
can be detected individually at SNR > 8. Some squares correspond to the 
same system at different times in its evolution. For illustration, the solid stars 
and circles represent the harmonics of a DWD with masses Mi = 0.74Mq 
and M2 = O.86M0 at two different stages of its evolution: at P = 872.64 s 
and e = 0.33 (stars) and P = 2211 .84 s and e = 0.74 (circles). In both cases 4 
harmonics ai'e observable individually at SNR > 8, but all n = 1 , . . . , 20 har- 
monics are plotted for illustration. The amplitudes are for a 5-year mission, a 
GC distance of lOkpc, and a single Michelson observable. 



TABLE 3 
Observable BINARIES. 



GC model 




Nohs at 3 kpc 


Nobs at 10 kpc 


Nobs at 50 kpc 


Typical 


29 


24 (28) 


17 (27) 


1 (7) 


47 Tuc-like 


102 


72 (96) 


36 (90) 


9 (24) 


Terzan 5-like 


140 


112(133) 


79(127) 


4 (36) 



Note. — Total number of eccentric LISA DWDs Mot and number of 
systems N^bs with at least two harmonics individually detectable by LISA at 
distances of 3, 10, and 50 kpc with a coherent SNR > 8 for a 5 year LISA 
mission. Both Not and N^bs are integrated over a 10-13 Gyr age range. Note 
that essentially all systems in GCs at Galactic distances are detectable with 
at least one harmonic; this number is provided within parenthesis in the Nabs 
columns. 

range from 24-112 for d = 3 kpc and 17-79 for d = 10 kpc. 
It is quite remarkable that LISA has a high enough sensitiv- 
ity to observe these systems throughout the Galaxy and, for 
a very limited region of the parameter space, even up to the 
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Large Magellanic Cloud. The probability to find at least 1 
detectable eccentric DWD is 1 % for the typical Galactic GC 
model and 3% for both the 47 Tuc-like and Terzan 5-like GC 
models (see Table |2]i. If we require only one harmonic to 
be detectable, the probabilities increase by a factor of ^ 2. 
Compared to the full population, the detectable samples of ec- 
centric DWDs furthermore favor higher orbital frequencies, 
total system masses, and chirp masses, and shorter system 
lifetimes; relevant orbital eccentricities span the entire range 
from e ~ 0. 1 to 0.8 (grey PDFs in Fig. [T]i. 

5. DISCUSSION 

In this Letter, we examine the formation of eccentric DWDs 
in the LISA band through dynamical interactions in GCs. Al- 
though the majority of cluster DWDs are circular, the exis- 
tence of this non-negligible population is in stark contrast to 
the Galactic disk population expected to consist exclusively 
of circular binaries. Owing to the non-zero orbital eccen- 
tricities, the DWDs emit GWs at multiple harmonics of the 
orbital frequency. For a 5-year LISA mission, we find a non- 
negligible number of systems for which multiple harmonics 
are detectable individually at SNR > 8, allowing the determi- 
nation of the orbital eccentricity from the ratio of the ampli- 
tudes of the detected harmonics. 

As a proof of principle, we show that for standard stellar 
and binary evolution assumptions, 47 Tuc and Terzan 5 type 
clusters are the most efficient eccentric DWD factories. For 
these two types of GCs, there is an 80% probability that at 
least one eccentric DWD is present in the LISA band at any 
given time. For a typical Galactic GC, this probability is 32%. 
Moreover, for a GC at lOkpc, the probability that one or more 
systems are present with at least two harmonics individually 
detectable at SNR > 8 is 1 % for a typical Galactic GC, and 
3% for 47 Tuc and Terzan 5 type GCs. Considering there 
are 79 known GCs within lOkpc from the Sun (Harris 1996), 
we thus expect at least a handful of eccentric DWDs to be 
detectable by LISA. There are furthermore 3 GCs within 2- 
3 kpc and 15 GCs within 5 kpc from the Sun. For these clus- 
ters, the SNR increases by a factor of ~ 2-3. These numbers 
may evidently vary significantly with the adopted stellar and 
binary evolution model parameters. Since our goal was to as- 
sess the feasibility of forming and detecting eccentric DWDs 
in GCs, we chose parameters which are common in standard 
population synthesis models. 

In principle, the independent knowledge of GC distances 



from electromagnetic observations allows the derivation of 
chirp masses of GC binaries from the amplitude of their 
GW signals. Unfortunately, the identification of sources 
with GCs is severely hampered by LISA'S limited angu- 
lar resolution. In particular, the error box in the sky is ~ 
1 .6(5 mHz/j/)2 (10/SNR)2 deg^ (e.g. Takahashi & Seto 2002) 
compared to a typical GCs size of < 0.002 deg^ Indepen- 
dent chirp mass and source distance determinations therefore 
rely completely on measurements of intrinsic changes of the 
GW frequency during the observation time (e.g. Takahashi 
& Seto 2002), assuming the drifts are entirely due to gravi- 
tational radiation with no significant contribution from tidal 
and/or magnetic spin-orbit coupling (e.g. Stroeer et al. 2005). 
It is therefore interesting to note that, depending on their phys- 
ical parameters and the phase of evolution at which they are 
observed, the systems in our simulations can produce observ- 
able frequency drifts. 

We conclude that the identification of eccentric binaries 
with LISA does not necessary imply association with NSs. 
Since binary component masses are difficult to derive from 
GW data, a clear discrimination between NS-WD or NS-NS 
and DWD binaries from LISA data may not be possible. How- 
ever, the number of eccentric NS-WD binaries in the three GC 
models is an order of magnitude lower than the number of 
eccentric DWDs, and NS-NS binaries are even scarcer Sta- 
tistically, an eccentricity detection therefore favors a DWD 
interpretation over a NS-WD or NS-NS interpretation. The 
formation of eccentric DWDs in GCs thus provides an ex- 
citing opportunity to study WD physics inaccessible through 
DWDs in the Galactic disk. Tidal effects particularly be- 
come increasingly important with increasing orbital eccen- 
tricity, either through the excitation of dynamical oscillations 
or through their effects on the orbital evolution. LISA could 
provide unique dynamical identifications of eccentric DWDs, 
opening up a unique astrophysical window to stellar and bi- 
nary evolution. 
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